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The Hierarchical Reference Theory (HRT) of fluids is a general framework for the description of 
phase transitions in microscopic models of classical and quantum statistical physics. The foundations 
of HRT are briefly reviewed in a self-consistent formulation which includes both the original sharp 
cut-off procedure and the smooth cut-off implementation, which has been recently investigated. The 
critical properties of HRT are summarized, together with the behavior of the theory at first order 
phase transitions. However, the emphasis of this presentation is on the close relationship between 
HRT and non perturbative renormalization group methods, as well as on recent generalizations of 
HRT to microscopic models of interest in soft matter and quantum many body physics. 



Liquid State Theory experienced remarkable successes in the understanding of the properties of dense, simple fluids 
[1] and more recently, of complex liquids, like colloidal suspensions or globular polymers 2]. However a fully quanti- 
tative description of phase transitions cannot be obtained by standard integral equations or diagrammatic expansions 
and requires specially devised approaches. The freezing transition, often dominated by geometrical constraints, is 
usually understood on the basis of Density Functional Theory [T], while fluctuation induced phenomena, like the 
liquid-vapour transition, need the completely different theoretical framework provided by the Renormalization Group 
method [3]. 

The Hierarchical Reference Theory of Fluids (HRT) was formulated more than 25 years ago as an attempt to 
reconcile the Renormalization Group (RG) approach, an extremely general and successful theory of second order 
phase transitions formulated in field-theoretical language, and the microscopic point of view adopted within Liquid 
State Theory [3] . The main goal of HRT was to develop a unified theory able to reproduce the universal properties at 
the critical point, as obtained by RG techniques, together with the non universal features of realistic models of fluids, 
present both in the critical region (e.g. critical density and temperature, short range correlations) as well as far from 
the critical point. 

A comprehensive presentation of the HRT formalism, including a detailed derivation of the basic equations and an 
extensive discussion of its critical properties, can be found in Ref. 5J. The earliest applications of HRT to simple 
fluids and to the Ising model are briefly discussed in the same reference. However, most of the more interesting results 
and extensions were developed later. By now it is therefore appropriate to summarize the advances in HRT which 
took place in the last decade, providing a guidance to the most significant recent applications. 

Such a critical discussion of the HRT formalism also gives the opportunity to emphasize its relationship with 
other theoretical tools which have been independently developed in different frameworks. Since the early nineties, 
sophisticated functional, non perturbative RG techniques (NPRG) were devised within field theory, leading to a boost 
of applications to many and diverse physical situations: from particle physics, to quantum many body theory, to 
classical statistical mechanics [S]- Only recently it was recognized that, despite the different language and physical 
context, a deep relation ties HRT and NPRG [7], whose developments have proceeded quite independently for 20 
years in distinct scientific communities. 



To introduce the Hierarchical Reference Theory of Fluids it is convenient to refer to a system of N classical particles 
confined in a volume V , in spatial dimension d, interacting via a pairwise potential v{v). Generalizations to classical 
spin systems, quantum fluids and magnets do not present additional difficulties and will be examined in the next 
Sections. The grand canonical partition function S of the model is written as: 



where the fugacity z is related to the chemical potential and to the thermal wavelength A = h/^/27rmkBT by 
z = A~'^ e^^ with f3 = {kBT)~^. An external field w(r) is also included for convenience and S is considered as a 
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functional of u{r). A detailed derivation of the HRT equations and a thorough discussion of their properties can be 
found in Ref. [5] via diagrammatic methods 0. Here we provide a brief introduction of the physical content of the 
method together with a simpler, less technical derivation. 

Since the pioneering works by Van der Waals, a large part of Liquid State Theory is based on the idea that 
the strongly repulsive short range interaction, distinctive of the interatomic potentials, provides a constraint to the 
available phase space, determining an important contribution to the entropy of the fluid, while the mainly attractive, 
regular part of i;(r) dominates the internal energy. This is certainly true for the special, although prototypical, case 
in which ^(r) has hard core at short distances, but holds also in the general case of soft core potentials. The different 
role played by the singular and regular contributions to the interparticle potential is manifest in lattice gas models, 
where particles live on the sites of a lattice: hard core repulsion between particles gives rise to the single occupancy 
constraint, which allows to map the lattice gas onto an Ising model with a spin-spin interaction proportional to 
the regular part of the interparticle potential. Mean field approaches developed for specific models, like Van der 
Waals equation in simple fluids or Weiss theory in magnets, suggest that only the latter interaction triggers the order 
parameter fluctuations, responsible for the order-disorder transitions. 

In fluid models, both the singular and the regular part contributions are included in the interaction v{r) which 
should then be split in a reference, singular, part u_R(r) and in the regular remainder w(r) = v(r) — vji{r). The partition 
function ([T]) can be written in a suggestive form by introducing the "configurational density" p{r) — '^^^{r — r^) 
whose thermal average is just the density profile p{r), which reduces to the ordinary average density p for u{r) — 0. 
In the following, the properties of the reference system will be considered as known. The fluctuating order parameter 
of the liquid-vapour phase transition in fluids is then identified as the configurational density and S is written as 

( -^ct>(0)\N r 

where the regular dimensionless interaction has been defined as 0(r) = —f3w{r) and the reference hamiltonian is 

= ^w/j(rj - Tj) + / dx-u(x)p(x) (3) 

According to the previous remarks, iffl, defines a non-trivial integration measure, limiting the phase space available 
to the particles. By factorizing the reference A^-particle canonical partition function Zpi = e~''"*«(^), Eq. ^ can be 
also expressed as: 

where the average {■ ■ ■) r refers to the Gibbs measure of the reference fluid, and the Fourier transforms of the potential 
0(q) and of the density fluctuations /5(q) — ^„ e*'* ''" have been defined. 

Such a form is particularly suitable as a starting point for a reformulation of Wilson's Renormalization Group 
picture in hamiltonian systems. In the momentum space RG formulation [3 [5], a sequence of systems, labeled by 
a wave- vector Q, is defined by allowing density fluctuations only at wave- vectors q larger than the cut-off Q. The 
evolution of the physical properties of the model is then followed starting from the limit Q = oo, where the system 
reduces to the reference one with a shift 'w{0)/2 in the chemical potential, to the final value Q — 0, where the cut-off 
is removed and the original partition function S is retrieved. Expression Q can easily accommodate the presence of 
a cut-off Q: 

EQ[u{r)]=J2 ^''''^p\ -^^^(^^ ^expQ J ^ ^(q) ^?(, - Q) p(q) p(_q)) (5) 

where 0{x) is the Heaviside step function. This definition embodies basic principle of HRT: A cut-off in density 
fluctuations is equivalent to a cut-off in the regular part of the two body interaction for the evaluation of the fluctuation 
corrections to the free energy, leading to the definition of a cut-off dependent potential WQ(q) ~ w{q) 9{q — Q). It is 
customary to perform a Legendre transform defining a cut-off dependent Helmholtz free energy Aq as a functional of 
the density profile p{r): 



Aq[p{v)] = -kT logSQ[7.(r)] - / dr (^(r) - kT logz) p{r) 



(6) 
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Finally, it is also convenient to include in the definition of the free energy functional the mean field contribution due 
to the residual interaction w{r) — WQ{r): 

AQ[pir)]^AQ[pir)] + ^ J drdr'[w{r-r')-WQ{v-r')]pir)p{r') (7) 

In this way, for Q — > cx), Aq reduces to the mean field approximation of the full free energy, which is recovered as 
Q — > 0. This way to artificially eliminate the long wave-length fluctuations from the thermodynamics is referred to as 
"sharp cut-ofF . However it is often more convenient to deal with a "smooth cut-ofF' procedure, where the 0-function 
in the definition of WQ{q) is substituted by a smoothed step function, denoted in the following as 6e{x) which is 
continuous and differentiable for every e > and tends to the step function in the e — > limit: 

WQ{q)^wiq)9,{q-Q) (8) 



A. The exact hierarchy in smooth and sharp cut-off formulations 

After having defined the cut-ofF dependent free energy, Aq via Eq. (l7| , the next step is to write the exact evolution 
(or flow) equation governing the change in the physical properties oithe model when the cut-off Q is varied. This 
goal is most easily achieved in the smooth cut-ofF formulation, where an infinitesimal increase of the cut-off Q leads 
to an infinitesimal change in the potential WQ{q). According to standard first order perturbation theory [1] we get, 
in a uniform fluid, 

15(-M,) _ p^ f^^^^^^^dj^^pjj^ 



V dQ 2 J ' dQ 2 dQ 

'^drF^(r)%M (10) 



2 J ' dQ 

where hQ{r) = gqi^) — 1 is the total correlation function of a fluid interacting through a potential WQ{r) and 
FQ{r) — pS{r) + p^/iQ(r). This smooth cut-off HRT equation has been obtained for the first time in Ref. [10] as an 
attempt to apply the Renormalization Group ideas to microscopic models. Note that the total correlation function 
hg rather than the radial distribution function gg appears in the evolution equation (j9| because the mean field 
contribution has been subtracted in the definition of the Q-dependent free energy ([t]). The total correlation function 
can be conveniently expressed in terms of the direct correlation function CQ(r) by use of the Ornstein-Zernike relation: 

The limit e — > in Eq. ([s]) requires some care, as already recognized in the description of the sharp cut-off HRT 
equation and also discussed in field theoretical approaches [TT] , because of the emergence of a discontinuity in hg (q) 
precisely on the same surface q = Q where the (5-function singularity develops in dgwg{q). Such a discontinuity 
originates from the Random Phase contribution to the exact direct correlation function, being the single term where 
the regular part of the potential appears outside any integration. If cj^(q) denotes the direct correlation function of 
the reference fiuid, the first terms of the perturbative expansion read 

cgiq) ^ cniq) + M^) + ■ ■ ■ (12) 

All the remaining contributions, being integrated, do not display discontinuities in the whole q-space. This suggests 
to introduce a modified correlation function Cg{q) which remains continuous also in the sharp cut-off limit e — )■ 0: 

CQ(q) = 5Q(q)-0Q(q) + 0(q) (13) 

Interestingly, this newly defined function is just the second functional derivative of the free energy functional 
—l3AQ[p{r)] with respect to density variations, evaluated at uniform density: 

J—, — TTT — V =CQ(ri-r2) (14) 

Sp{ri)Sp{r2) p(r)=p ' p 

where the last term is due to the ideal gas contribution, included in our definition of Ag. We first write the right 
hand side of Eq. (|9| by introducing the Fourier transforms: 

1 d{-f3Ag) _p f dq 1 d4>g{q) ^^^^ 



V dQ 2 y (27r)'i i„^CQ(q) + p 0(q)-0Q(q) 



dQ 
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The easiest way to take the sharp cut-off hmit of this flow equation is to introduce the differential operator Oq acting 
only on the Q-dependence of 0Q(q) (while CQ(q) is treated as Q-independent) : 



1 d{-(3AQ) 
V dQ 



J (0^^^^°^ [l-P (CQ(q) + <^Q(q)-0(q) 



Taking the e — > limit and writing dg as a finite difference, Eq. (16) is then written in the form [1]: 



1 dj-^AQ) 
V dQ 



dl7q 

(27r)' 



log 



1+J-Q(q),^(q) 



(16) 



(17) 



where the integration is limited to the spherical surface (in dimension d) defined by q — Q and the function Tq{c\) is 
defined by the Ornstein-Zernike relation as: 



p 



l-pCQ(q) 



(18) 



Equations ( 15 ) and ( 17 ) govern the HRT flow of the free energy in the smooth and sharp formulations respectively. As 
previously mentioned, the initial condition aX Q = oo coincides in both cases with the known mean field approximation, 
while for Q ^ all density fluctuations are included and the exact free energy is in principle recovered if ^Q(q) is 
known. 

HRT then provides the exact evolution of the free energy in terms of the two body correlation function Fq (q) (or 
equivalently J-q{c\)) whose flow is determined by similar equations involving the three and four particle correlations 
[31[inj- These equations, which can be derived along the same lines previously sketched, are most easily obtained by 
functional differentiation of ( jTsj ) and p?] ) with respect to p{v) leading to a hierarchy of differential equations describing 
the RG flow of the free energy and of the n-particle direct correlation functions c5^(ri • • • r„), formally defined as the 
^th functional derivative of {—PAq) with respect to the density profile p(r) evaluated at uniform density p(r) — p. The 
ideal gas contribution is included in our definition of Aq leading to an additional term {n — 2) \ (— p)^~" Oi^^i ^(^i — ri) 
with respect to the standard definition of the many particle direct correlation functions 0. 

As an example, we explicitly report the flow equations for the previously defined two point function in smooth 
cut-off formulation 10 



g5Q(k) 
dQ 



d(f'Qiq) 
dQ 



c?(k, q, -k - q) FQik + q) {-k, -q, k + q) F^iq] 



d4'Qi<i) 
dQ 



(19) 



and in sharp cut-off formulation [4]: 

dCgik) _ 1 /■ da 



^c^(k,-k, q,-q) 



•^|(q)'^(q) 

1+J-Q(q)0(q) 



dQ 2 J^^Q (27r)° 

/ ^ C?(k, q, -k - q) F,(k + q) c?(-k, -q, k + q)-^|^^ 
Jq^Q (Ztt)" l+J^Q(q)0(q) 

where, in the latter equation, FQ{k + q) is defined as: 



F (k + q) = J + [l+-^o(k + q)0(k + q)J for |k + q| < Q 

J"Q(k + q) for|k + q|>Q 



(20) 



(21) 



Note that, in the limit k — ^ 0, these flow equations reproduce exactly the second density derivative of the corresponding 
free energy evolution ( 15|17 ), showing that the density dependence of the free energy flow encodes some (partial) 
information contained in the full infinite hierarchy. 



B. HRT as a microscopic implementation of the Functional RG 



Since the first derivation of the HRT equations in the sharp cut-off formulation [5] a close similarity to some 
previous differential RG formulation was noticed. In particular it was shown that close to a critical point, in the long 
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wave-length limit the full HRT hierarchy precisely reproduces the flow equations for the one particle irreducible RG 
generator [12]. The analogy between the Hierarchical Reference Theory and functional RG approaches is however 
even deeper, as recently shown by Caillol [?]• To highlight the connection between the two approaches it is convenient 
to start from the definition of the partition function Sg (|5|). As pointed out in an early attempt to apply RG 
techniques to microscopic models [13 , by performing a Hubbard- Stratonovich transformation, the partition function 
can be expressed as a functional integral over an auxiliary scalar field ipq governed by an effective action Sq[iPci]- 
This commonly adopted procedure, however, is rather formal and does not lend a well defined physical meaning for 
the Hubbard Stratonovich field V'q- It is possible to circumvent this problem by first introducing a reference effective 
action SR[ipq\ by: 



log 



N ' '' q 



(22) 



This definition clearly forces the effective field V'q to faithfully represent density fluctuations. The exact partition 
function is then written as 



Sq^ = ^ E'^Q(q)V'qV'^q+ ^"^^"^ 



V ^\ 2 
q 



+ /3u(q) V-q 



(23) 
(24) 



This expression shows that S'gfV'q] can be written as the sum of the effective action in the absence of any cut-off S'['0q] 
(i.e. the value attained at Q = 0) plus a "momentum dependent mass term" [Bj: 



(25) 



with 



^Q(q) 



q) 



(q) [l-0,iq-Q)] 



(26) 



Following Wetterich we now perform a "modified Legendre transform" H] introducing the field pq conjugated to 
the external potential u{(l), which, due to the definition (24) coincides with the statistical average of the field ipq- 



r[pq] = -logSQ -f E(^(^) " P-^ ^ ^ E^Q(^)^q^-q 

q 



(27) 



Comparing this definition to Eq. ([?]) it is clear that the average action coincides with the modified free energy defined 
in HRT: F = 13 Aq, providing the link between HRT, based on the gradual inclusion of the Fourier components of the 
regular part of the physical interaction, and the "effective average action" RG method built upon the idea of adding 
a cut-off and momentum dependent mass term Rq (q) to the original action. From Eq. ( 26 1 it follows that for Q = 
the mass term vanishes identically, while at fixed Q, Rq{ci) is large only ioi q < Q proviHing an effective cut-off for 
the long wave-length fluctuations, as requested in all the implementations of the momentum space RG ideas [9 . A 
difference in the initial condition however persists between the average action method and HRT: in the RG approach, 
Rq{(1) diverges as Q ~> oo, thus providing an infinite mass to all fiuctuations and leading to a soluble problem. 
Instead in HRT we get limQ_j.oo i?Q (q) = 0(q): for this special but finite value of the mass term the regular part 
of the interaction cancels exactly in Sqlipq] and the model reduces to the soluble reference system. The procedure 
adopted in HRT is reminiscent of some recent development in functional RG methods [13]. 



C. Approximate closures 

The solution of the infinite hierarchy of equations previously derived needs some approximation scheme. Two kinds 
of approximations are customary in this type of problems: truncations and closures. A truncation of the hierarchy 
simply corresponds to neglecting n-particle correlations for n larger than a maximum value n* (typically n* =4). 
Such a procedure has been followed in the early investigations of RG equations but it is rigorously justified only near 
a critical point in spatial dimension close to d = 4 [3], while it is not expected to provide an accurate description 
of a system, like a liquid, where short range correlations are important. Instead, most of the applications of HRT, 
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both in the smooth and sharp cut-ofF formulations, adopted a non perturbative closure to the first equation of the 
hierarchy ( 15|17 ) and only this type of approximation will be discussed here. Non perturbative closures in fact allow to 
preserve the convexity of the free energy, providing a physically satisfactory description of phase transitions between 
homogeneous phases, together with a picture of the critical region fully consistent with the scaling hypothesis. 

The evolution equation for the free energy involves the knowledge of the two particle direct correlation function in 
the presence of a cut-ofF Q. A closure is a parametrization of the direct correlation function expressing Cq{\C} in terms 
of the free energy density Aq /V . As previously noticed, the compressibility relation, which follows directly from Eq. 
pl. 



V 



CQ(k-O) 



1 

p 



(28) 



plays a special role in HRT because, when implemented into the free energy evolution equation, allows to include in 
the very first equation some information contained into the whole hierarchy. 

The simplest non perturbative closure of the first HRT equation is inspired by the well known Random Phase 
Approximation (RPA) in Liquid State Theory [T]: 



CQ(k)=C^(k)+AQ0(k) 



(29) 



The exact initial condition at Q — >■ oo is obtained for Aq — 1, while the compressibility relation (28 1 can be enforced 
at all Q via a suitable choice of the parameter Xq. When this closure is substituted into Eq. (15) or (17) the resulting 
partial differential equation for the free energy Aq as a function of density p and cut-ofF Q can be solved for a given 
value of the temperature. 



Interestingly, this approximate closure can be equivalently formulated starting from the evolution equations ( 15 ) 



or (171 for the full functional Aq[p{v)] [4]: By restricting its form through a suitably chosen parametrization, the 
functional flow equation gives rise to partial differential equations for the parameters. For instance, the RPA-like 
closure is obtained by use of a simple Weighted Density Approximation 15j for the free energy: 



Aq[p{v)]^Ar[p{v)]+ / dr/Q(p(r)) 



(30) 

where ^^[r] is the exact free energy functional of the reference system and the function fqip) is the parameter whose 
evolution will be determined by the flow equation. The weighted density in Eq. (301 p(r) is defined by the local 
average: 



p(r) = J dr'a;(r — r') p(r') 



(31) 



with a normalized weight function given by uj{q) = [0(q)/0(O)]^/^ [16 . By substituting this form into the evolution 



( |l7| , we obtain the same expression found by use of the ansatz (29) if we identify Aq w{0) 



equation ( 15 ) or 

This class of closures, in the NPRG framework, are known as Local Potential Approximation (LPA) because of the 



_ d'fQ 



ap2 



similarity between Eq. ( 30 ) and the usual Local Density Approximation of Density Functional Theory. 



Clearly, more refined approximations to the direct correlation function, or equivalently to the free energy functional, 
may be devised, leading to more accurate representation of the exact free energy flow. In fluids with a hard sphere 
reference system (in the continuum or on a lattice) , Eq. ( 29 ) can be generalized by including a further contribution, 
non vanishing only inside the hard sphere diameter, so as to satisfy the core condition gg^r) = 0. This constraint can 
be implemented exactly in the smooth cut-ofF Formulation for Yukawa potentials [17] thanks to the exact analytical 
solution of the Ornstein-Zernike equation for this model |18) . 

Analogous closures might be implemented on the flow equation for the two particle direct correlation function 
(19) or (p0|. A closure to the second equation of the hierarchy amounts to expressing the three and four particle 



correlations in terms of the two point function. A possible closure has been proposed already in the earliest stage of 
HRT [H [TU] but it has not been implemented yet: 

, _ aCQ(k + q) 



c?(k, q. 



(k, -k,q,-q) = 



dp 

dp^ 



g^CQ(k-q) 
dp^ 



(32) 



(33) 



This class of approximations has the advantage of allowing for a non trivial anomalous dimension exponent at the 
critical point [H [10] . A similar closure has been independently proposed in the context of functional RG [TS' and a 
simplified form has been numerically investigated in the same framework |20j . 



Up to now, all the applications of HRT to microscopic models adopted a closure of the form ( 29 ) with or without 



the additional optimization of the potential inside the hard sphere diameter to enforce the core condition. 
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D. HRT in Smooth cut-ofF and the Self Consistent Ornstein Zernike Approximation 

The Self consistent Ornstein Zernike Approximation (SCOZA) is a well known liquid state theory which proved 
successful in the determination of phase diagrams in fluids characterized by a pair interaction written as the sum of 
a singular hard sphere part plus a Yukawa tail w{r) [21j. Originally, SCOZA has been formulated as a generalization 
of the Mean Spherical Approximation [T] for the direct correlation function outside the core: 

c(r) = A0(r) (34) 

while enforcing the core condition g{r) — to define c(r) inside the hard sphere diameter. The parameter A is 
determined by requiring thermodynamic consistency between the internal energy and the compressibility routes to 
thermodynamics: 



dp' 



2 



'h-i 



,2 

dr g{r) w{r) 



d_ 
'dp 

92 f-BAn\ r , , . 1 



drc{r)-- (35) 



The two exact expressions (35), together with the core condition and the parametrization (34) outside the hard sphere 



diameter, give rise to a partial differential equation for the free energy as a function of p and (3. Remarkably, this 



equation is exactly reproduced by HRT in the smooth cut-off formulation ( 10 ) by use of a linear turning on of the 
interaction 

Wt{r) = tw{r) (36) 

where t is the switching-on parameter [T7|. Clearly this form does not satisfy the requirements stated in Section [Tl| 
for a selective introduction of density fluctuations: while SCOZA just tunes the amplitude of w(r), HRT changes 
simultaneously both amplitude and range of the potential by imposing a cut-off Q on its long wave-length Fourier 
components thereby qualifying as a microscopic implementation of RG ideas. In fact, SCOZA has been shown to 
provide a good description of the non universal properties in the critical region, as well as non classical critical 
exponents but fails in reproducing the correct scaling picture predicted by RG and HRT . An attempt to reconcile 
HRT and SCOZA has been formulated in Ref. [33], where a further adjustable parameter is added in the RPA-like 
closure of the HRT equation, to enforce thermodynamic consistency between the compressibility and the internal 
energy route. 

III. CRITICAL PROPERTIES AND PHASE COEXISTENCE 

One of the goals which motivated the formulation of HRT was to provide a description of the critical properties 
of realistic fluid models consistent with the RG picture. Critical exponents, governing the singularities of thermo- 
dynamic quantities in the critical region, acquire universal values quoted in Table |l] [24j . At the critical density, the 
constant volume specific heat diverges as Cy t~°' as the reduced temperature t — T/T^ — 1 vanishes; the isothermal 
compressibility as K^t"^ and the correlation length as ^^t~'^. The shape of the coexistence curve is asymptotically 
given by \p{T) — pc\ \t\^ , while the equation of state on the critical isotherm is given by \P — Pc\ -^\p — Pc\^ . Finally, 
precisely at the critical point, the pair correlation displays a power law tail of the form h{r) ~ 

In the sharp cut-off formulation, universality is indeed a direct consequence of the asymptotic form of the HRT 
hierarchy close to the critical point |3]: all the microscopic details of the model disappear from the asymptotic 
equations and the full hierarchy acquires a universal form equivalent to a previously studied RG differential generator 
[12j . Scaling laws and non classical values for the critical exponents immediately follow from this identification. In 
the smooth cut-off formulation, the universality of critical phenomena is less apparent because the specific form of 
the chosen cut-off function WQ(r) ( 8|26 ) survives also in the asymptotic critical region. Universality of the critical 



properties can then be proved only order by order in the e = 4 — d expansion and no general proof is available 1101 
Here we will briefly review the critical properties of the approximate non perturbative class of closures previously 



introduced both in the sharp and in the smooth cut-off formulation ( 29 ) . Moreover we will sketch the argument 
leading to the conclusion that already at the level of this simple approximation, HRT always gives a convex free 
energy, even at coexistence. The unphysical divergence of the compressibility at the phase boundary predicted by 
the sharp cut-ofF closure |25| is circumvented by the smooth cut-off formulation which then provides a satisfactory 
description of first order transitions . 
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Exponent 
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= C+/C- 


"Exact" 


0.110 


0.327 


1.237 


0.630 


4.789 


0.036 


4.76 


HRT-sharp 


-0.07 


0.345 


1.378 


0.689 


5 








HRT-smooth 


0.05 


0.330 


1.300 


0.650 


5 
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TABLE I: HRT estimates of tiie critical exponents and compressibility amplitude ratio in three dimensions for the sharp and the 
smooth formulations of HRT compared to the exact values [24] obtained by extrapolation of high-temperature series expansions. 
The smooth cut-off procedure adopted here is detailed in Section |III B| 



Other, more sophisticated closures have been proposed. As previously noticed, the second equation of the HRT 
hierarchy describes the evolution of the two particle direct correlation function in terms of the three and four particle 
correlations ( 19|20 1. By use of the closure ( 32|33 ) it has been shown that the critical exponents are correctly reproduced 
to second order in e = 4 — d, including a non vanishing anomalous dimension ?7 [11 [5l [10] . Unfortunately this quite 
promising non perturbative closure has not been thoroughly investigated in microscopic models, probably due to the 
problems posed by the numerical solution of the resulting non linear partial differential equation. Recently, a similar 
(although simpler) approximation has been studied in the asymptotic regime, in the framework of NPRG, leading to 
encouraging results for the universal quantities |20j . 



A. Sharp cut-off 



In the sharp cut-off HRT, the free energy evolution equation ( 17 1 depends on the two point function J-q (q) evaluated 
"on shell", i.e. on the surface q = Q. This means that in thelate stage of the integration, when Q is small, only 
the long wave-length physics contributes to the free energy flow. Moreover, in the critical region, J-Q(q) diverges at 
small q and then we can disregard the additive term in Eq. |l7[ ) and the sub-leading log 0(q) contribution. Finally, 
the evolution equation acquires a universal form 



1 g(~Mg) 
V dQ 



(37) 



where Kd is a geometric factor equal to the ratio between the volume of the unit sphere in d dimensions and {2tt)'^. 
When we insert a closure of the form (291 at right hand side, recalling Eq. (18 1, we need to evaluate the direct 
correlation function at wave-vector Q ^ 0. By exploiting the analyticity of c/j(q) and (/)(q) around q ~ 0, we can 
write 



l-pCQ{0) + pb( 



(38) 



where b determines the curvature of the direct correlation function in this approximation and attains a finite, positive 
limit at the critical point. Recalling the compressibility sum rule (28 1 and inserting this asymptotic form into Eq. 



(37) we obtain a closed partial differential equation expressing the free energy flow at long wave-length and in the 
critical region. The RG structure of this equation can be highlighted by a suitable rescaling, as shown in Refs. [HIS], 
which allows the explicit calculation of the universal properties, like critical exponents and scaling functions, via the 
fixed point analysis. The analyticity in the momentum dependence of the direct correlation function, implied by the 
RPA-like ansatz (29), immediately gives the vanishing of the anomalous dimension, or Fisher exponent, 77 = in any 



spatial dimension. This result is exact for d > 4, when the critical exponents are classical, while for d < 4 such a 
simple non perturbative approximation to the first HRT equation gives critical properties correct to first order in the 
e = 4 — d expansion. However this closure leads to a qualitatively incorrect prediction in d = 2, when the anomalous 
dimension is crucial in ensuring the existence of a fixed point solution. Few critical exponents in d = 3 are reported 
in Table |T] together with the accepted values obtained from extrapolations of the high temperature expansion in the 
Ising model [24) . 

It has been realized in Ref. [25] that the class of non-perturbative closures to the first HRT equation previously 



introduced (29) gives rise to a convex free energy in the limit Q — 7> 0: long wave-length fiuctuations are responsible 



for the flattening of the isotherms in the coexistence region. A similar conclusion has been reached in the framework 
of NPRG in Ref. [27 . An analytical argument justifying the convexity of the free energy can be provided starting 
from the asymptotic equation ( 37 ) previously derived. Below the critical temperature, mean held approximation gives 
negative compressibility in a finite region around the critical density p^, and mean field approximation acts as initial 
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condition to the HRT flow equation, which describes the effects of fluctuations in reshaping the free energy. Eq. ( 38 ) 
shows that J-q{Q) grows large at small Q. In this regime it is convenient to introduce the auxiliary quantity 



uqip) = - log Tq{Q) = log 



V 



which satisfies the flow equation 



dQ 



(39) 



(40) 



When J-q{Q) is large, uq becomes large and negative and the exponentially vanishing term at left hand side of Eq. 
( 40 1 can be neglected. Then Eq. ( 40 ) shows that the only possible asymptotic form of uq (p) is 



2b 

die,. 



[ip - p,f ~ A^] Q 



2-d 



(41) 



where A is an integration constant. However, this solution is consistent with the assumption uq{p) — t- —oo only for 
Pc — A < p < pc + A (and d > 2): below the critical temperature, the compressibility diverges within a density interval 
(coexistence region) of width A. Clearly, this asymptotic solution breaks down on the binodal, i.e. for p — pc = ±A. 
As discussed in Ref. [2S], a study of Eq. (40) near the binodal shows that just outside the coexistence curve the 
inverse compressibility behaves as 



V dp^ 



cx 



\P- Pc 



A)a^ 



(42) 



for d < 4 leading to an unphysical divergence of the compressibility at the phase boundary. Interestingly, above four 
dimension the correct behavior is restored and the compressibility attains a finite limit on the coexistence curve. 



B. Smooth cut-ofF 



The analysis of the critical properties within the smooth cut-off formulation of HRT is slightly more subtle than 



in the sharp cut-off case. Formally, the free energy evolution equation ( 15 ) depends on the two point function at 



all wavelengths, contrary to the sharp cut-off case, and some preliminarily consideration on the cut-off depend ence 
of the function 0Q(k;) is appropriate in order to simplify the asymptotic flow. According to the identification (26) 
the difference i?Q(k) — (/)(k) — (?!)Q(k) acts as a "mass term" in the associated field theory: as such, it should scale 
as times a short range function of the ratio k/Q 6 . In previous applications of HRT to microscopic models 
[iniEHlISS] the choice i?Q(k) = e~^* 4){ke^) was adopted at small values of Q: the momentum cut-off is proportional 
to e"'* with s : Q ^ oo, while the cut-off function is identified as the physical two body interaction. In this way an 
efficient suppression of long wave-length density fluctuations is achieved, together with the correct asymptotic scaling. 
Moreover, this special choice allowed to set the initial condition of the HRT flow at s = 0, when '?iQ(k) identically 
vanishes and the physical properties of the model reduce to those of the known reference system. Instead, when HRT 
was applied to the study of the 0^ scalar field theory ;26J the simpler form 



i?Q(k) = ^(k) - 0Q(k) - 5 (0^ _ fc2) ^(g _ ^) 



(43) 



was adopted. As usual b defines the long wave-length behaviour of the direct correlation function Cq (k) ^ Cq (0) — bk^. 
Such a form of the cut-off function Rq (k) , closely resembles the choice proposed in Ref. [30 and allows to highlight 
the main features of the simple RPA closure to HRT in the smooth cut-off formulation, avoiding many of the technical 
inconveniences of other options. By substituting Eq. (43) into the HRT flow equation (15) we obtain [26] : 



1 dj-PAg) 
V dQ 



= -bKdQ 



d+l 



dp^ 



The initial condition is set at the ultraviolet cut-off of the 0"* theory, taken as the inverse length unit {Q — 1): 

fiAi[p) _ pAiipc) 



V 



V 



+ t{p- Pc) +U{p- Pc) 



(44) 



(45) 
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FIG. 



1: Inverse compressibility (48 1 of the 4>* theory in d = 3 below the critical temperature. The dashed line corresponds to the 
initial condition, which coincides with the mean field approximation, while the solid line is the result at convergence (Q = 0). 
Data refer to u = 0.035 and t = —0.25. Upper (lower) panel shows the results of the sharp (smooth) cut-off formulation of 
HRT. 



Here i is a measure of the mean field reduced temperature t — {T — Tmf) /Tmf while u is the self- interaction defining 
the bare theory. It is convenient to rescale the free energy by Kd and the density by \f~KdJh to get rid of the 
dimensional constants in the evolution equation (44). A further rescaling of both the density and the free energy 



by suitable powers of Q confers a RG structure to this flow equation, whose fixed point corresponds to the critical 
singularity. A standard stability analysis of the fixed point solution allows to evaluate the critical exponents, whose 
values in d = 3 are shown in Table |l] Similarly to the sharp cut-off case, the parametrization ( |29[ ) provides the correct 
exponents to first order in the e = 4 — c? expansion. Being Cq (k) analytic in the wave- vector, the anomalous dimension 
vanishes (77 = 0) for all d as in the sharp cut-off case. 

In Refs. [3T] it was suggested that the smooth cut-off formulation of RG qualitatively improves on the description 
of the system close to the phase boundary. This expectation was later confirmed both within HRT [35] and NPRG 
[32j . Following a procedure similar to that adopted in the sharp cut-off case, we first define the auxiliary quantity 



(46) 

By writing the flow 



which, below the critical temperature, diverges to -l-oo as Q — > within the coexistence region 
equation for uq and equating the diverging terms to leading order we get the solution 

UQ{p)={^^-{p-p,f\Q-'' (47) 

where Pc and A are integration constants. This solution is consistent with the assumption uq — > -t-00 for \p — pd < A 
showing that pc and A are the critical density and the width of the coexistence curve respectively. Therefore, also the 
smooth cut-off formulation of HRT gives rise to a convex free energy. Moreover, a detailed analysis of the solution in 
an infinitesimal neighbourhood of the phase boundary p ~ pc ± A shows that a discontinuity develops for Q — >■ [55] . 
This behaviour is clearly illustrated in Fig. [l] which shows the evolution of the (running) inverse compressibility 

(48) 



contrasting the sharp and the smooth cut-off flows. As suggested in Ref. 



the presence of a discontinuity at the 

The 



phase boundary is a generic feature of any smooth cut-off scheme and not a peculiarity of the chosen form ( 43 1 
finite value attained by the isothermal compressibility x the phase boundary allows to introduce the compressibility 
amplitude ratio defined as U2 ~ Cj^ /C_ where the ± label refers to the case above or below the critical temperature. 
If T is the reduced temperature t = (T — Tc)/Tc, 



x{p*)-c±{±Ty 



(49) 



For r > the density is fixed at its critical value p* = pc while for r < the compressibility is measured along the 
phase boundary: p* ~ Pc + A(t). This amplitude ratio U2 is known to be universal and the smooth cut-off HRT 
estimate agrees well with high temperature extrapolations, as shown in Table |T] The full equation of state in scaling 
form can be also evaluated as shown in Ref. [23 ■ 
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FIG. 2: Snapshots of the inverse compressibiUty as a function of density at three values of the cut-off Q for u = 0.035. Left 
panel refers to t = —0.25 below tc ~ —0.203989, right panel to t = —0.20 above the critical temperature. Dotted line (red): 
Q — 1 (which coincides with mean field approximation); dashed line (blue): Q = 0.135; full line (black): Q = 0. In the right 
panel the effects of fluctuations for Q < 0.135 are negligible making the dashed and full lines almost indistinguishable. 



C. Locating the spinodal within HRT 



The spinodal is a useful way to discriminate between intrinsically unstable and metastable states, even if this curve 
is not rigorously defined within equilibrium statistical mechanics, the free energy being a convex function of density 
in the whole phase diagram. In Ref. [29] it was noticed that a special feature of the smooth cut-off HRT flow below 
the critical temperature suggests a way to discriminate between instability and metastability. At mean field level, 
a negative compressibility identifies the unstable region lying inside the spinodal curve. However, fluctuations are 
expected to modify the instability boundary which indeed evolves as Q lowers, as shown in Fig. [2] short wave-length 
fluctuations up to a characteristic length-scale Rc slightly shrink the instability region, while fluctuations at wave- 
vector Q ^ Qc = lead to a widening of the two-phase region, defined by the asymptotic vanishing of the inverse 
compressibility. Therefore the HRT flow allows to identify a density interval, deep inside the coexistence region, 
characterized by a negative compressibility at all wave-lengths: the unstable region bounded by the spinodal line. 
Instead, close to the phase boundary, but still inside the binodal, a region exists where the compressibility remains 
positive even by including fluctuations up to a maximum length-scale Rc' the metastable region. The curves dividing 
the regions of positive and negative Xq^ i'^ the (p, Q) plane at different temperatures are shown in Fig. jsj together 
with the extent of the metastable region obtained via the criterion previously outlined. The cut-off value Qc where 
the width of the unstable region has a minimum, vanishes when the critical point is approached. This corresponds 
to the divergence of a characteristic length-scale Rc ^/Qc, which may be identified as the critical droplet radius on 
the spinodal line, according to the standard droplet model picture of nucleation [53] ■ This divergence follows a power 
law consistent with the scaling Rc^ ^ ^ I''' I" with = 7/2 (recall that within the present approximation the critical 
exponent 77 vanishes). 



IV. APPLICATIONS TO MICROSCOPIC MODELS 



RPA-like closures of the type discussed in Section p^I Cj have been applied to several microscopic models both in the 
continuum and on the lattice within the sharp cut-off formulation of HRT. The earliest studies were focused on simple 
liquids, like noble gases or Lennard- Jones fluids [31] . The RPA closure ( 29 ) was later suitably generalized to include the 
correct short range behavior of correlations characterizing the liquid phase, the so called core condition, which forces 
the vanishing of the radial distribution function inside the hard sphere radius. Extensions of this procedure to deal 
with soft core potentials, were also devised [5]. The Ising model on a cubic lattice was investigated as the inescapable 
test case of any theory which qualifies for the description of phase transitions [35] • More recent applications, together 
with the first implementations of the smooth cut-off formulation of HRT and an extension to binary fluids are briefly 
reviewed in the following Sections. 



12 




FIG. 3: Boundaries of the region defined by a negative compressibility as a function of the cut-off Q for u = 0.035. At Q = 1 
the standard mean field approximation is recovered, while all fluctuations are included at the top of the figure. Left panel 
refers to t = —0.25 < t^. At each Q the shaded area is the metastable region, where the compressibility becomes negative due 
to the inclusion of fluctuations at wavelengths Rc < \ < , where _Rc ~ 3 at this temperature. Dots mark the unstable 
region, characterized by a negative compressibility at all wavelengths A < . At convergence (top of the figure) full dots 
mark the coexistence region (located at \p — pc\ ~ 1.5), empty dots the "spinodal" (at \p — pc\ ~ 1.05). Right panel, same plot 
for t = —0.20 > tc- Here the unstable region closes at Qc ~ 2.45, and the fully interacting system is in a single phase state. 



A. One component fluids 



The very first physical system studied by HRT (in the sharp cut-ofF formulation) was a microscopic model for simple 
fluids [34^ , which allowed for a detailed comparison with experiments in rare gases both for the thermodynamics and 
for structural data [3Bj. Later, HRT was also applied to complex fluids, a model for Cgo |3Z|i and more recently to 
systems of interest in soft matter physics. An HRT code in smooth cut-off formulation has been also developed and 
applied to the study of a Hard Core Yukawa fluid for different values of the inverse range z [17] and to the Screened 
Restricted Primitive model of electrolytes [31] . Here we briefly review three physical systems where the HRT in sharp 
cut-off formulation has been successfully applied: a colloidal suspension with added polymers, a system of colloidal 
particles with competitive interactions and a model of star polymers in solution. The aim of this presentation is 
to show the flexibility of the HRT formalism in dealing with markedly different classes of systems. The set up of a 
numerical code for the solution of the HRT non-linear partial differential equation presents some technical problem 
due to the occurrence of singularities below the critical temperature, related to the divergence of the isothermal 
compressibility in the whole coexistence region. These problems, already discussed in the literature [ST] |39l |40| , are 
also briefly illustrated in the Appendix. 

Depletion interactions are known to be ubiquitous in soft matter physics. The most celebrated system where the 
depletion mechanism is able to induce particle aggregation is a suspension of colloids (modeled as hard spheres) where 
some non-ionic depletion agent, usually a polymer, is added. The depletant is treated as an ideal gas, while an 
excluded volume interaction acts between polymers and colloids. The effective colloid-colloid pair interaction turns 
out to be attractive and its analytical form can be calculated exactly in this model [2J: 



3 r 



1 



2{l + q) 2(1 + 9)3 



for Kr <l + q (50) 



Distances are measured in units of the colloid diameter a and q = 2Rg/a, where Rg is the gyration radius of 
the polymer. This effective interaction is proportional to the (reservoir) polymer packing fraction rjp and, being of 
entropic origin, to the temperature. HRT in sharp cut-off formulation has been tested against numerical simulations 



for this model [41] using the hard sphere gas as reference system and the Asakura-Oosawa (AO) potential (50) as 



attractive tail. The very short ranged potentials which occur quite often in soft matter physics cannot be accurately 



represented by use of the RPA-like non-perturbative closure (29 1 adopted in most of the implementations of HRT. 
Therefore, in the study of the Asakura-Oosawa model moderately large size ratios, q > 0.4, have been considered. 
For these systems many body interactions, which appear for q > 0.154, become relevant and affect the location of 
the coexistence curve, as shown in the numerical study performed in Ref. [41 . However, the comparison between 
simulations and HRT results for the AO model shows a good agreement and allows to quantify the extent of the 
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FIG. 4: HRT results in sharp cut-ofT formulation for the AO model (50 1. Upper panel: true correlation length ^ as a function 
of the reduced temperature t = -S^ — 1 in a logarithmic plot. Its slope defines the critical exponent v. Lower panel: effective 



critical exponent Ueff, normalized to its asymptotic HRT value u 
Core- Yukawa fluid (full squares) . 



0.69, for the AO model (open circles) and for a Hard 



critical region, also characterizing the crossover between a mean-field classical pre-critical region and the asymptotic 
regime. A nice power-law divergence of the correlation length in an extended range of reduced temperatures is shown 
in the upper panel of Fig. |4j while in the lower panel the temperature dependence of the effective critical exponent 
is compared to the case of a standard Hard Core- Yukawa fluid with z — 1.8, which is known to mimic the properties 
of simple fluids. Interestingly, the AO model displays a power-law behaviour in a significatively larger domain of 
reduced temperatures [42 . 

Charged colloids repel each other via a screened Yukawa potential which, together with an attractive short ranged 
attraction, due either to dispersion forces or to the presence of a depletion agent in solution, give rise to an effective, 
tunable, potential characterized by interparticle attraction at short distances and a repulsive tail at longer range. 
Few experimental results [15], together with a theoretical analysis based on simulations [H], suggest the possible 
occurrence of modulated phases of different morphologies. A class of potentials with these features was investigated 
via the sharp cut-off HRT |,45j in the homogeneous fluid phase, where the incipient cluster formation may induce 
peculiar thermodynamic and structural anomalies. The pair interaction of the model is of the form named Hard Core 
Two Yukawa Fluid (HCTYF): 

v{r) = \ f ^ J (51) 

— e \- Ae tor r > 1 

\ r r 

As usual, distances are measured in units of the hard core diameter. The reference system is defined by the hard 
sphere part plus the repulsive Yukawa tail, while w{r) is identified with the attractive contribution. The investigation 
in Ref. |45j has been carried out for the choice zi = 1 and Z2 = 0.5. Temperatures are measured in units of the energy 
scale e while the parameter A > is varied. Within HRT the fluid phase is stable against microphase formation 
for A < 0.097 but even such a small repulsive contribution to the interparticle interaction gives rise to remarkable 
changes to the phase behavior of the fluid. When the parameter A approaches the stability limit, the coexistence 
curve considerably flattens near the critical point and a large region of enhanced compressibility opens above the 
critical temperature. The same picture is confirmed by use of SCOZA, which provides a phase diagram topology in 
quantitative agreement with HRT. In Fig. [5] SCOZA results are shown to display such a remarkable flattening of 
the binodal. The physical region for such a peculiar behavior may be identifled as the presence of large clusters of 
strongly correlated particles, which anticipate the formation of modulated phases. This can be seen in Fig. |6] where, 
in the upper panel, the presence of two correlation lengths below a threshold temperature is illustrated. While a single 
correlation length does diverge at the critical point, which belongs to the Ising universality class, the second, "intra- 
cluster" , correlation length remains remarkably large in the whole critical region. The presence of two correlation 
lengths reflects in the occurrence of anomalies of the effective critical exponent on approaching the critical point, as 
shown in the lower panel of Fig. [6]for A = 0.097. A sharp growth of both correlations above the critical temperature 
gives rise to strong peaks and sharp jumps in Veff, while the true asymptotic behavior of the correlation lengths can 
be seen only very close to the critical temperature. 

Star polymers in a good solvent display an effective interaction characterized by an ultra-soft repulsive core of 
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FIG. 5: Coexistence regions of the HCTYM, obtained by SCOZA, for two choices of the amplitudes of the repulsive tail: 
A = (red line) and A = 0.097 (blue line). The dashed lines limit the regions of enhanced compressibility, defined as the 
thermodynamic states where the dimensionless isothermal compressibility S{0) is lager than 10. 
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FIG. 6: HRT results in sharp cut-off formulation for the HCTY model |45] at ^ = 0.097. Upper panel: true correlation length 
^ as a function of the reduced temperature t = — 1 in a logarithmic plot. The diverging correlation is shown by full circles. 
Open squares represent the "intra-cluster" correlation length. Lower panel: effective critical exponent Veff, normalized to its 
HRT value f — 0.69, for the asymptotically diverging correlation. 



entropic origin, which can be conveniently parametrized as: 



VRir) 



Inr 



fcT A/3/2 (,^^ 



2 



-1 1 
e 2 



v7(i-'-) 



for r < 1 
for r > 1 



(52) 



where distances are measured in units of the corona diameter of the star a and the functionaUty / of a star polymer 
coincides with the number of its arms. This expression accurately reproduces Molecular Dynamics simulations for 
/ > 10 f^H. Here we will discuss the case / = 32, displaying a stable fluid phase at all concentrations p* (in units of 
(T^^). HRT has been used to study the consequences of an additional attractive interaction between star polymers, on 



top of the entropic repulsive contribution (52): such an attraction may be induced by different physical mechanisms, 
like the presence of residual dispersion forces or the effects of a depletant added in the solution. Recent simulations 
P7] confirmed that the attractive interaction between star polymers due to added depletant can be conveniently 
parametrized by a Fermi function; 



wir) = -C 



1 



(53) 



15 




12 3 4 

P 



FIG. 7: HRT coexistence curve for the model of star polymers in solution discussed in the text [IS]. A non-linear temperature 
scale has been used to enhance the visibility of the binodal structure (the temperature scale has been set for convenience to 
Tq = 0.23). Note the presence of two critical points: one at extremely low concentration (p*i = 0.074, T*i = 0.6123) and the 
other at hquid densities (p*2 = 1.952, r*2 = 0.2471). The triple point is located at (p* = 0.975, T^* = 0.2426). 



This model was previously studied by means of the sharp cut-off formulation of HRT [IB] for / = 32 and the choice 
of the parameters A = 2.1, B — 0.35, while C played the role of a temperature scale which defines the dimensionless 
temperature T* = kT/C. For this model the reference system, defined by the potential (52) 
and its thermodynamic and structural properties, which serve as an input to the HRT forma : 



is highly non trivial 
, ism, were obtained by 

use of accurate integral equations of Liquid State Theory. In the investigation carried out in Ref. [15] the Modified 
Hypernetted Chain equation ;I was adopted. The resulting phase diagram is shown in Fig. Tjand displays two critical 
points separated by a triple point. Both the extremely low critical concentration p*^ and the very presence of the 
high density phase transition are consequences of the ultra-soft nature of the repulsive effective interaction. 



B. Lattice models 



Applying HRT to lattice models requires a few preliminary modifications of the evolution equations due to the 
different nature of the order parameter and to the geometry of the Brillouin Zone (BZ). Let us consider an interacting 
spin model on a lattice, defined by the hamiltonian 

if = - ^ J(R-R')Sr-Sr, -/i^^^ (54) 

R,R' R 

where are n-componcnt unit vectors defined on the sites R of a lattice. The spin coupling J(R) and the external 
magnetic field h in the z direction are the parameters defining the model. We will consider the case of un-frustrated 
ferromagnetic systems (J(R) > 0), where the order parameter of the transition can be identified as the magnetization 
per site and can be written as a thermodynamic derivative of the partition function Zih): 

^ ^^aiogz(M 

- ^^^^ - N dph ^^^> 

In anti-ferromagnetic models ( J(R) < 0) it would be appropriate to introduce a symmetry breaking staggered external 
field instead of h. According to the basic principles of HRT, we define a sequence of interactions Jq(R), labeled by an 
infra-red cut-off Q which suppresses the long wavelength Fourier components, thereby inhibiting the order parameter 
fluctuations on those scales. As usual, it is convenient to introduce the cut-off dependent free energy Aq^ui) by 
performing a Legendre transform on log Z{h). Then we define a modified free energy Aq by including the full mean 
field contribution ([7|: 

Aq (m) = (m) - ^ ^ [ J(R - R') - Jq (R - R')] (56) 

R,R' 
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The free energy evolution equation at fixed temperature and magnetization m (along the z axis) can be immediately 
obtained by first order perturbation theory: 



1 dAq 
N dQ 



R 



dJQ{^) 
dQ 



(57) 



where -F'q(R) is the connected spin-spin correlation function 

FQ(R-R') = (SR-SR,>Q-m2 (58) 
and the canonical average {■ • •)q is taken at fixed magnetization m and interaction Jq(R). Eq. (58) can be written 



in terms of the longitudinal i^Q(R) and transverse Fg (R) correlations 

i^^(R-R') = {SI,S^,)q-' 
Fq(R-R') = {S^S^,)q 



with the result 



1 OAq 

N dQ 



dk 



4(k) + (n-l)F^(k) 



dQ 



IBZ (2^)" 

where the integral is restricted to the appropriate Brillouin Zone of the chosen d-dimensional lattice. 



(59) 
(60) 

(61) 
(62) 



Following the same steps detailed in Section 



II A 



first define new correlation functions by including in i<g"(k) 



the mean field contribution due to the remaining part of the interaction J(k) — jQ(k): 

F""(k) 



J-r(k) 



l-/3F^"(k) J(k)-jQ(k) 



(63) 



Then we recall the thermodynamic relations expressing the longitudinal and transverse susceptibilities in terms of the 
corresponding correlation functions at vanishing momentum: 



= /3 



92 (Aqim) 



dm? 



N 



N 



And finally, we parametrize these correlation functions in a simple RPA-like form, analogous to Eq. ( 29 ) 

hT 



(64) 



(65) 



A" - J(k) 

which defines a non-perturbative closure to the evolution equation, provided the parameters A" and A"*" are determined 
by the exact sum rules (64). Analogously to the previously examined continuum case, it is convenient to define the 

where the 0-function selects a surface 



(66) 



cut-off dependent interaction JQ(k) by an expression of the form of Eq. 
where the interaction J(k) is constant. In the sharp cut-off limit, a possible choice is: 



('j(k) for J(k) < J(0) (1 - Q) 
1 elsewhere 



For (3 = 1 — Min J(k)/J(0) the system reduces to non-interacting spins. By lowering Q, more and more Fourier 
components of J(k) are included into the free energy until, in the Q — >■ limit, the full interaction is recovered. With 
such a choice, the HRT flow equation for the free energy (57) in sharp cut-off acquires the form: 



1 dAq 
N dQ 



kT 



D[Q) 



log 



f /3./(0) 

g/3j(o) 



+ {n — 1) log 



[X 



-Ll-l 



g/3j(o) 



(67) 
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Exponent 


n = 


n = 1 


n = 2 


n = 3 


HRT-sharp 


0.606 


0.689 


0.768 


0.826 


"Exact" 


0.588 


0.630 


0.67 


0.71 



TABLE II: The correlation length critical exponent u predicted by Eq. ( |67| l on the basis of a fixed point analysis of the HRT 
equation in sharp cut-off formulation. The "exact" values are theoretical estimates [2j: obtained by combining extrapolations 
of high-temperature series expansions and field theoretical resummations. The limiting value u — 1 for n — >■ cxd is correctly 
predicted by HRT. 



where D{Q) is the "density of states" defined by 

dk 



The critical properties provided by this equation are exact to first order both in the e = A — d and in the 1/n expansion. 
The HRT results for the correlation length exponent v in three dimensions, obtained by numerically solving the fixed 
point equation associated to Eq. (67), are shovifn in Table [TT| for few values of n. In this approximation the two point 



correlations are analytic functions of and then the anomalous dimension exponent rj vanishes, while the theoretical 
estimates suggest r] ~ 0.03 — 0.04 for all the values of n considered here. The other HRT exponents can be obtained 
by use of the scaling and hyperscaUng relations which are rigorously verified within this approximation. Notice that 
in the sharp cut-off HRT formalism, the critical exponents directly follow from the single assumption of analvticitv 
of [J^Q(k)] ^ on k^ at all Q's 5/. any approximation retaining this property is bound to the values shown in Table |ll| 
which constantly overestimate by about 10% the accepted values. HRT is shown to provide the best results for n = 0, 
where the n-vector model describes the self avoiding walk problem, while their quality somehow worsens at larger n, 
although they reproduce the correct limit v — 1 as n — ;> cxd. 

The Ising model on a cubic lattice was the first lattice system investigated by sharp cut-off HRT [3S] by use 
of a more sophisticated closure, belonging to the class previously discussed, which consistently preserves the spin 
normalization constraint = 1. A thorough analysis of lattice models via the smooth cut-off version of HRT 
has not been attempted yet, although a very similar procedure, put forward in the framework of non perturbative 
Renormalization Group, has been recently developed |14j . The results shown in |35j provided a very accurate estimate 
of the critical temperature of the Ising model and compared favourably with the coexistence curve predicted by high 
temperature expansion. A numerical study of the same HRT equation aimed at determining the effective parameters 
of a scalar field theory able to reproduce the correct pre-critical behaviour of the Ising model was performed in Ref. 
[33] ■ Interestingly, HRT showed that a (j)* representation of the long wave-length spin fiuctuations was inadequate for 
a faithful representation of the non-universal properties of the Ising model, while the inclusion of a (j>^ term allowed 
to considerably improve the description of the critical region of the model. 

In Ref. [5U] the HRT formalism was generalized to the case of an uniaxial spin model where can assume three 
values: —1, 0, 1. This system can be easily mapped into a symmetric binary mixture on a lattice with opposite signs 
for like and unlike interactions. If the coupling is limited to nearest neighbours, the Blume-Capel model is obtained, 
while for long range interactions decaying as ~ 1/-R, the Lattice Restricted Primitive Model (LRPM) is recovered. 
The phase diagrams of both models in the symmetry plane (zero magnetization) were obtained by integrating the 
HRT flow equation [50 . The topologies of the phase boundaries in the density-temperature plane, shown in Fig. [sjfor 
the LRPM, turn out to be very similar: at high temperatures a critical A-line separates a charge-ordered fluid phase 
at low density from a homogeneous plasma phase. The critical line ends at a tricritical point where it merges into 
the liquid-vapour coexistence curve. The comparison with available simulation data for the LRPM shows that the 
structure of the phase diagram and the universality classes of the transitions are correctly reproduced, although some 
quantitative discrepancy is present, suggesting that a better representation of pair correlations should be devised. 



C. Binary Mixtures 
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The HRT formalism can be easily generalized to the case of binary mixtures |53j : the flow equations in the smooth 
and sharp cut-off (17 1 formulation remain valid for mixtures provided the two body interaction 0(q) and the 



pair correlations ^b(q) and Cq^q) are substituted by matrices whose indices run over the two species. The Ornstein- 
Zcrnike relation ( 18 1 acquires the form 



rife 



(69) 
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FIG. 8: Phase diagram of the LRPM in the (p, T) plane. HRT results: the tricritical point is marked by an open star. Full 
dots give the A-line above the tricritical temperature and the liquid-vapour binodal at lower temperatures. The SCOZA data 
for the A-line are shown by the dotted line. Simulation data from Refs. [511 [52] are represented by filled triangles and open 
squares respectively. 



while the compressibility sum rule (28) becomes a set of three independent equations vifhich can be written as: 

92 



dpidpj 
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Pi 



(70) 



The derivation of the HRT evolution equations straightforwardly follows the same steps sketched in Section jll A| By 
adopting the smooth cut-off procedure, the resulting expression is: 



1 d{-PAQ) 
V dQ 
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where matrix products are understood at right hand side and "Tr" stands for the trace operation. In the sharp cut-off 
limit this form reduces to 
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(72) 



( 72 ) must be interpreted in matrix sense and coincides with the 
q)0(q) . By use of a RPA-like closure of the form (29 1, 

the 



The trace of the logarithm at right hand side of Eq. 
logarithm of the determinant of the 2x2 matrix ! 

the flow equation then becomes a partial differential equation in three independent variables: the cut-off Q and 
densities of the two species (pi, ^2)- 

The analysis of criticality in binary mixtures is performed following the procedure developed in simple fluids j5j 153) : 
several universality classes are predicted by the HRT equations, but most of the possible fixed points turn out to 
be unstable and the critical properties at a generic point along a transition line is of the Ising type. Actually, two 
distinct fixed points correspond to critical exponents related to the case of a one component scalar order parameter. 
One is weakly unstable and gives rise to strictly Ising- like critical properties, while the other is stable and correctly 
reproduces the so called phenomenological approach by Fisher which, starting from the concept of field mixing, 
predicts the occurrence of some weak exponent renormalization |54| : the long wavelength limit of a structure factor 
Sij{Q) at the critical density and concentration diverges with an exponent 7/ (1 — a) (with a ~ 0.11), a 12% larger than 
the corresponding value for one component fluids, while the isothermal compressibility displays a weak divergence 
governed by the exponent a /{I — a). HRT allows to identify the thermodynamic perturbation corresponding to the 
crossover between these two fixed points |55j . thereby providing a microscopic interpretation of the experimental 
difficulty in observing the true asymptotic critical behaviour. Special points in the phase diagram might display 
genuinely different critical indices, like un-renormalized Ising or tricritical exponents. 
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HRT also naturally suggests an unbiased definition of the order parameter at a critical point of binary fluids 
|53j . which in general cannot be simply identified either as density or concentration, because critical fluctuations are 
expected to affect both quantities. The critical point in binary fluids is signaled by the divergence of the leading 
eigenvalue of the (symmetric) matrix of the structure factors Sy^k) for k — ^ 0. The structure factors are a direct 
measure of density fluctuations and can be expressed as the equilibrium average of a product of density operators 
Sij(k) = jj (/5i(k)/5j(— k)). To discriminate between demixing and liquid- vapour transitions in binary mixtures it 
is convenient to introduce suitable linear combinations of SijCk), the Bathia- Thornton structure factors |56j . which 
properly characterize the density and concentration fluctuations: 

Spp = 5n(k) + 522(k)+2 5i2(k) 

Sec = S'ii(k) + 522(k) - 2x1 X2 S'i2(k) 

Spc = 5'cp = a;i S'22(k) - a;2 S'ii(k) + (xi - X2) S'i2(k) (73) 

where Xi — pi/ p, p — pi+p2 are the number concentrations of the two species and the total density, respectively. Let us 
denote by A leading eigenvalue of the 2x2 matrix built with the previously defined Bathia Thornton structure factors 
in the k — limit, and (cos 6*, sin 6*) be the corresponding normalized eigenvector. According to these definitions, 
the diverging eigenvalue provides a measure of the fluctuations of the particular linear combination of density and 
concentration ip(k) = cos0p(k) + p sin6'a;2(k) (for k — )■ 0): limk-j-o ('0(k)'0(— k)) = A, suggesting the identification 
of ip as the order parameter of the transition. It is convenient to pictorially represent the order parameter as an arrow 
in the {p, X2) plane pointing in the direction identified by 9. Remarkably, the "mixing angle" 9 is related to well 
defined thermodynamic quantities [S]: 

1 - a;2 cot^ = pvi (74) 
^ ptaii9 . (75) 



A(l/p) 

The first equation is written in terms of the partial molar volumes v^ = (dV/dNi)p^T which attain a finite limit at 
the critical point. The second, equivalent, relation expresses 9 in terms of the differences of concentration Ax2 and 
molar volumes A{l/p) at coexistence on an isotherm-isobar near the critical point. Clearly, the direction of the order 
parameter is defined at each point of a critical line and varies continuously along that line. While in some mixture, 
the character of the transition remains almost unaltered by varying the thermodynamic state, in other cases it may 
instead experience considerable changes, passing from demixing to liquid-vapour by reducing the density. A plot of 
the mean field critical lines for two Lennard- Jones binary fluids mimicking two rare gas mixtures, shown in Fig. |9j 
illustrates these two possibilities. For convenience, the total packing fraction 7/ — f (picr^ -|- P20'2) and the volume 
fraction 772/'? are used in place of total density p and number concentration X2. 

The numerical integration of the HRT equation poses some challenge due to the presence of three independent 
variables: cut-off Q, density p and concentration X2- Special attention must be paid to the stability requirements and 
to the possible occurrence of singularities induced by the presence of large coexistence regions in the phase diagram. 
A stable numerical code has been developed for two classes of models of binary fluids: the rare gas mixtures ^57^ and 
the symmetric Yukawa mixtures [58] , which include the Ising fluids, defined as hard spheres with an internal Ising-like 
internal degree of freedom whose relative orientation specifies the sign of the pair interaction |59j . In both cases, the 
HRT equation in the sharp cut-off formulation has been numerically integrated. Here we just review few results for 
the model of rare gas mixtures, where extensive experimental data are available [50" ■ 

The Argon-Krypton system is modeled as a mixture of particles interacting via Lennard- Jones (LJ) potentials 
[57j . Each LJ interaction has been preliminarily mapped into a reference hard sphere potential plus a LJ attractive 
tail. The effective hard sphere diameter does not coincide with the LJ parameter a and is determined according to 
the standard Barker-Henderson prescription [I] which gives rise to a weak temperature dependence. The bare LJ 
parameters adopted in the numerical calculation have been fixed to reproduce the location of the critical points of 
the two pure components: CTAr = 3.52 A, ctki = 3.77 A, eAr-Ar = 113.8 K and CKr-Kr = 157.9 K. The strength of the 
unlike attractive interaction has been set to eAi— Kr = 133.5 K, and additivity of the effective hard sphere diameters 
has been assumed. Plots of the coexistence regions in the packing fraction/molar concentration plane are given for 
two temperatures in Figs. |10||1H where the tie-lines, connecting thermodynamic states in equilibrium, are also shown. 
The Argon-Krypton mixture belongs to the topology of phase diagram illustrated in the left panel of Fig[9j where the 
critical points of the two species are connected by a line of critical points. A different view of the coexistence regions 



at the same temperatures can be seen in the molar concentration/pressure plane in Fig. 12 at the higher temperature 
pure Argon is in the fluid phase, while at the lower temperature the coexistence regions of the two components merge 
and define a unique, connected domain. A good agreement between theory and experiments j61j can be appreciated 
for this system. Finally, it is instructive to check that, within HRT, the full coexistence region collapses to a single 
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FIG. 9: Mean field critical lines projected in the plane total packing fraction/Krypton volume fraction, for two parameter 
choices mimicking an Argon-Krypton (left panel) and a Nenon-Krypton mixture (right panel). The arrows identify the mixing 
angle 9 defined in the text: an arrow parallel to the density axis indicates a liquid-vapour transition, an arrow along the volume 
fraction axis, a demixing transition. 




FIG. 10: Coexistence regions in the packing fraction/Krypton molar fraction in the Argon-Krypton mixture according to HRT 
at T = — 95.77°C. Tie lines in the two-phase regions are shown. In the coexistence regions darker colors correspond to higher 
pressures. 



line when intensive thermodynamic quantities are selected, as required by thermodynamics: in Fig. [13] the same 
coexistence regions are plotted in a diagram where the difference between the chemical potentials of the two species 
and the pressure are considered as independent variables. 

The phase diagram of the Neon-Krypton mixture instead displays a different topology, belonging to the class 
depicted in the right panel of Fig. [9] The same Lennard- Jones model used for the Argon-Krypton mixture has 
been adopted with parameters appropriate for the Neon-Neon (ctng = 2.85 A, eNe-No = 33.5 K) and Neon-Krypton 
interactions (eNe-Kr = 54.4 K). The coexistence regions of the two pure components remain separated in the whole 
phase diagram: while pure Neon is in the gas phase in a large temperature range, the Neon-Krypton mixture displays 
phase coexistence between states of considerably different composition. The order parameter smoothly changes 
character signaling a liquid-vapour transition at low Neon concentration and a demixing transition at higher pressures 
as already captured by the mean field calculations of Fig [9j The coexistence region in the three planes packing 



fraction/molar concentration, molar concentration/pressure and chemical potential/pressure are shown in Figs. 14 



15 16 and |17| respectively. The comparison with the experimental data [62 shows good agreement at not too high 



pressure, while some discrepancy is present at high pressures, where the repulsive part of the interaction plays an 
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FIG. 11: Coexistence regions in the packing fraction/Krypton molar fraction in the Argon-Krypton mixture according to HRT 
at T = — 130°C. Tie hnes in the two-phase regions are shown. In the coexistence regions darker colors correspond to higher 
pressures. 
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FIG. 12: Coexistence regions in the Krypton molar fraction/pressure plane for the Argon-Krypton mixture according to HRT. 
Left panel: T = — 95.77°C, right panel T = — 130°C. Experimental phase boundaries [61] are shown as red dots. 

important role. As discussed in Ref . [37] , it is likely that a better representation of the soft core part of the interatomic 
potential must be employed for a fully quantitative description of demixing transitions at high pressures. 

V. HRT FOR QUANTUM SYSTEMS 

HRT was originally developed in the framework of classical statistical physics but its formalism can be generalized 
to quantum hamiltonians. Lattice spin models have been considered in Ref. [53], where the reference system was iden- 
tified as free spins in an external magnetic field which couples with the order parameter. Quantum antiferromagnets 
on an hypercubic lattice in d dimensions are described by the isotropic, nearest neighbour Heisenberg hamiltonian 

H = J J2 SR-SR,-/i,^e^s-^5^ (76) 



and require the introduction of a staggered external field hs- Here, Sr, are spin S operators defined on the sites R 
of the lattice and g = (tt, tt . . . ) is the antiferromagnetic vector in d dimensions. The reference system is defined by 
the hamiltonian ( 76 ) for J — and gives rise to a solvable, single site problem. Following the standard sharp cut-off 
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FIG. 13: Coexistence regions for the Argon-Krypton mixture according to HRT plotted in the plane dimensionless chemical 
potential difference/pressure. Left panel: T = — 95.77°C, right panel T = — 130°C. Note the collapse of the full coexistence 
region into a single line. 
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FIG. 14: Coexistence regions in the packing fraction/Krypton molar fraction in the Neon-Krypton mixture according to HRT 
at T = — 110°C. Tie lines in the two-phase regions are shown. In the coexistence regions darker colors correspond to higher 
pressures. 



formulation of HRT, the spin-spin interaction is gradually turned on, by including its Fourier components starting 
from the boundary of the antiferromagnetic Brillouin Zone. The exact evolution equation for the free energy per site 
in sharp cut-off formulation was obtained in Ref. [63j and involves the two body dynamical spin-spin correlations in 
imaginary time: 



1 OAq 
N dQ 



kT 



E 



dk 



J(0)2 - J(k)2 - q\ { log [l - J--(k,C.„) J(k) 



2 ^Jszi^nr 

log [(l - J-^-(k, c.„) J(k)) (l - J-3-(k + g, c.„) J(k + g)) + J-^^(k, c.„) J^^y{k + g, J(k) J(k + g 



(77) 



where J(k) is the Fourier transform of the spin-spin nearest negighbour interaction and Q ranges from J(0), where 
fluctuations are absent and the mean field approximation is recovered, to Q = 0, where all Fourier components 
are introduced. The structure of this equation is rather similar to its classical counterpart (67) with two main 
differences: the presence of off diagonal two body correlations J^Q^(k,a;„), due to the non commutativity of the 
spin components, and the summation over the Matsubara frequencies w„ = 2t: kTn. The imaginary time variable 
in fact naturally emerges when quantum perturbation theory is employed to evaluate the effects of an additional 
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FIG. 15: Coexistence regions in the packing fraction/Krypton molar fraction in the Neon-Krypton mixture according to HRT 
at T = — 107°C. Tie lines in the two-phase regions are shown. In the coexistence regions darker colors correspond to higher 
pressures. 
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FIG. 16: Coexistence regions in the Krypton molar fraction/pressure plane for the Neon-Krypton mixture according to HRT. 
Left panel: T = "110°C, right panel T = — 107°C. Experimental phase boundaries [HI] are shown as red dots. 



infinitesimal momentum shell in the spin-spin interaction. Notice that in this formulation, the cut-off on fluctuations 
is set only on the spatial Fourier components, while the Matsubara frequency variable runs on the whole real axis 



for all Q's. The Quantum HRT (QHRT) equation appropriate for the antiferromagnetic Heisenberg model ( 76 ) was 
studied by employing a simple closure similar in spirit to the classical RPA-like approximation ( |65[ ): the structure of 
the dynamical correlations was evaluated at mean field level and only "mass renormalization" was allowed by imposing 
thermodynamic consistency between the zero frequency limit of the dynamical correlations at wave-vector k = g and 
the derivatives of the free energy Aq with respect to the staggered magnetization. The explicit form of the adopted 
closure reads: 



- ^(k) 

TO7^a;2 + /4 - J(k)2 



= -2 , . . (78) 



J-o^(k,c^) = 5 \ (79) 

-^Q (k,c.) = -^rhr, (80) 
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FIG. 17: Coexistence regions for the Neon-Krypton mixture according to HRT plotted in the plane dimensionless chemical 
potential difference/pressure. Left panel: T — — 107°C, right panel T — — 110°C 
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FIG. 18: QHRT results for the antiferromagnetic Heisenberg model in d = 3 for spin S — oo and 5* — 1/2. Left panel: QHRT 
spontaneous magnetization for 5 = 1/2 and S = oo (classical limit). Right panel: constant field specific heat at ft = as 
a function of temperature compared with classical (full dots) [B?] and quantum (full squares) ;65: Monte Carlo simulations, 
respectively. In the = oo limit the magnetization is normalized to S and the specific heat to . 



where is the staggered magnetization, and the two cut-off dependent parameters (/ij^,/i|| 
transverse and longitudinal staggered susceptibilities: 



d (Aq 
dms V N 



= rus n±_ + J(g) 



are related to the 

(81) 
(82) 



In Fig. [18] the constant field specific heat and the spontaneous magnetization of a Heisenberg antiferromagnet on 
the cubic lattice below the critical temperature are shown for two choices of the spin S: the extreme quantum value 
5 = 1/2 and the classical limit S* — >■ oo. Monte Carlo simulations for Ch [MlISS] are also reported. A more detailed 
comparison between theory and simulations can be found in ^53!, showing that the performances of QHRT have 
generally a lower quality with respect to the previously shown applications to classical models, suggesting that a 
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better representation of static and dynamic spin correlations has to be developed: the different role of longitudinal 
and transverse correlations is probably not fully grasped by the adopted single mode approximation, which did not 
prove accurate enough to provide a consistent picture of the crossover between the renormalized classical and the 
quantum critical regime of the Non Linear Sigma Model fBS] . 

The effects of easy-axis anisotropy have been also investigated by use of QHRT 67j , as well as the role of a weak next 
nearest neighbour interaction [68]. However, a QHRT analysis of quantum (i.e. zero temperature) phase transitions, 
or of the occurrence of off-diagonal long range order has not been attempted yet. As shown in the related non 
perturbative RG approach for the Bose-Hubbard model [55] , the investigation of Bose condensation by RG methods 
requires a different cut-off procedure and it may be more convenient to set a cut-off on the kinetic energy term rather 
than on the two body interaction. According to this scheme, the reference system becomes a purely classical model, 
while quantum fluctuations are gradually introduced, leading to the appropriate free energy flow equation. 

The smooth cut-off formulation of HRT has not been applied to quantum systems yet, although the derivation of the 
flow equations does not present any additional difficulty. A closer connection between QHRT and Non Perturbative 
Renormalization Group methods may provide some clue on the structure of more refined closures to the flow equation. 

VI. CONCLUSIONS 

A proper description of first and second order phase transitions within liquid state theory requires a careful account 
of long wavelength fluctuations. The Hierarchical Reference Theory of fiuids was developed precisely to achieve this 
goal. While the general framework of the theory was set since the very first formulation of HRT, the implementations 
to specific physical systems required the development of approximation schemes which may vary according to the 
individual features of the model. In this review we tried to summarize the efforts in this direction carried out in the 
last ten years. Among the features shared by all the non-perturbative closures of the HRT equation we may mention 
the correct scaling form of the thermodynamics and correlations in the critical region, the prediction of non-classical 
critical exponents and the exact implementation of the Maxwell construction at first order transitions, showing that 
the inclusion of long wave-length fiuctuations gives rise to coexistence regions at constant pressure and chemical 
potential. The latter feature is particularly valuable in the case of mixtures, where the ad hoc implementation of 
the equal tangent construction, required by most liquid state theories, may prove rather cumbersome. The smooth 
cut-off formulation of HRT, proposed several years ago, has been quantitatively investigated only recently. Among its 
peculiar properties, it is worth mentioning the absence of the spurious divergence of the isothermal compressibility at 
the phase boundary: a problem which affects the sharp cut-off formulation. This allows the location of the spinodal 
loci within a Renormalization Group approach, as explictly shown for the Yukawa fluid. Extensions of the smooth 
cut-off HRT equations to the case of mixtures or quantum systems are straightforward and represent a fertile research 
fleld for future studies. 

The approximations to the HRT hierarchy we have been considering so far were limited to closures to the flrst 
equation, which provides the flow equation for the free energy of the model. More sophisticated closures, at the level 
of the second equation, which contains information on the two particle correlations, have been formulated but not 
quantitatively investigated yet. 

Soft matter physics provides a stimulating environment for any microscopic theory able to capture the mechanisms 
leading to phase transitions. And HRT was applied to several systems of physical interest in this context. However, 
we feel that improved analytical representations of pair correlations in different frameworks have still to be developed 
to achieve a degree of accuracy comparable to simulation methods. For instance, a common feature of colloidal 
suspensions is the presence of short range attractive interactions, often represented by the "sticky limit" . A closure of 
the form ( |29| ), adopted in most of the implementations of HRT, implies a linear relation between the direct correlation 
function and the attractive part of the potential, which is known to provide poor results for short range interactions 
[T]. Extensions of HRT to lattice systems, mixtures and quantum models have been briefly covered in this review, 
but other interesting directions were explored in the literature, as the application of HRT to polymer physics or to 
disordered media [70]. 

A special mention must be devoted to the relationship between HRT and non perturbative (or functional) Renor- 
malization Group methods. The two approaches are in fact intimately related, as already noticed by Caillol [7], 
and have been independently developed in different physical frameworks. We hope that the present review, which 
speciflcally devoted Section [HB] to this problem, will contribute to clarify this issue providing a dictionary between 
the two, substantially equivalent, formulations. 

It is a pleasure to express our gratitude to all the former students and collaborators who helped us in these 
investigations, most notably D. Pini, M. Tau and F. Lo Verso. We acknowledge financial support by the CIMAINA 
foundation (LR) and by MiUR through a PRIN/08 grant (AP). 
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VII. APPENDIX 



The flow equation for the free energy (17), together with the definition (18) and the parametrization of the direct 



correlation function (29), supplemented by the compressibility sum rule (28 1 form a closed partial differential equation 
(PDE) for the free energy as a function of the density p and the cut-off Q. As a practical example here we report 
the explicit equations for a model defined by a spherically symmetric interaction written as the sum of a reference 
potential vji{r) (and corresponding free energy Ar{p)) and an attractive tail w{r). As usual we set (j){r) = —/3w{r)] 
Fourier transforms are denoted by tildes. 

The initial condition of the PDE is set at some large value of Qoo, such that the pair interaction w(q) can be 
neglected for q > Qoo ■ The free energy at Qoo is then given by the mean field approximation: 



(83) 



Here and in the following, the ideal gas contribution is included in the definition of Aj^. The integration of the 
PDE is then performed backwards, down to Q = 0, where all fluctuations are included and Aq coincides with the 
free energy of the model. A mesh is defined on the density axis {0,pmax)- Here Pmax is some relatively high value 
(typically Pmax ~ 1 in units of the hard-sphere diameter) and the free energy Aq (p) is propagated from Q to Q — dQ 
by discretizing the PDE on the mesh. It is convenient to write the PDE ( [l7| in quasi-linear form by introducing the 
auxiliary quantity 



uq{p)^ log l + TQ{Q)4>iQ) 



The evolution equation for this new variable can be straightforwardly obtained from Eq. (|T7| 
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where as usual Kd = {6Tr^)~^ in three dimensions (d = 3) and cii(Q) is the direct correlation function of the reference 
system. Usually the hard sphere gas is adopted as a reference system and c^XQ) is evaluated in Percus-Yevick [T] or 
Verlet-Weis approximation [7T]. This PDE has the structure of a non-linear diffusion equation, where —Q plays the 
role of time-like variable and p of space-like variable. The stability requirement for this class of PDE's is given by the 
positiveness of the effective diffusion coefficient. This condition is verified by Eq. (85) provided (f){Q) is positive. 

Boundary conditions at p = and p = Pmax a-re required: it follows from the definition that at lo w d ensity 
uq{p) ~^ p4'{Q)i while UQ^pmax) is usually approximated by its mean field expression, i.e. by Eq. (84) with 
!Fq{Q) = p/[l — PCr{Q) — p4>{Q)]. Different boundary conditions at high density give very similar solutions to the 
PDE, at least when the range of the potential is not too short. 

Details on the discretization procedure can be found in Refs. [331 [SJ. We finally remark that at low values of Q 
(e.g. Q < Qo ^ 0.1) it is convenient to adopt a logarithmic scale in the cut-off, by setting Q = Qo e^* with t £ (0, oo) 
to slow down the inclusion of long wave-length fiuctuations. 
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